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A fully gauge- invariant (pseudoscalar) meson photoproduction amplitude off a nucleon including 
the final-state interaction is derived. The approach based on a comprehensive field-theoretical 
formalism developed earlier by one of the authors replaces certain dynamical features of the full 
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allows for a systematic improvement of this approximation. The feasibility of the approach is 
illustrated by applying it to both the neutral and charged pion photoproductions. 
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I. INTRODUCTION 

Ever since the pioneering work by Chew, Goldberger, 
Low, and Nambu on pion production [1], the study of 
photo- and electroproduction of mesons off nucleons has 
been utilized as one of the major research avenues to 
learn about the excited states of the nucleon. In order to 
extract accurate information on nucleon resonances, one 
needs — in addition to precise and extensive experimen- 
tal data — reliable reaction theories that allow one to 
disentangle the resonance contributions from the back- 
ground contributions to the observables. 

The extant descriptions of meson photoproduction 
reactions span a wide range of different approaches 
(e.g., Chiral Perturbation Theory, tree-level effective 
Lagrangians, K- matrix approach, etc. [2-10]). The 
present work is based on the dynamical framework of 
meson-exchange models of hadronic interactions [11-14] 
in which the composite nature of hadronic vertices is ac- 
counted for by so-called form factors. Since, at present, 
our theoretical understanding of these vertex form fac- 
tors is rather incomplete, one usually parameterizes the 
vertex structure by phenomenological functions (usually 
of monopole or dipole form) whose parameters are ad- 
justed to fit the data. The presence of such form fac- 
tors spoils gauge invariance of the photoproduction am- 
plitude already when dressing the bare tree level in this 
phenomenological manner. The inclusion of an explicit 
hadronic final-state interaction (FSI) further complicates 
this problem since the construction of the corresponding 
interaction current (where the photon interacts with the 
hadronic structure within the vertex) consistent with the 
FSI requires the knowledge of the underlying dynamical 
structure, a requirement that is impossible to be satisfied 
in an approach where phenomenological form factors are 
employed. In order to maintain gauge invariance in this 
situation, one needs to resort to finding prescriptions that 
are consistent at a phenomenological level with the vari- 
ous dynamical models. The existing prescriptions [15-22] 



do not, and indeed cannot, provide a unique answer to 
this problem since manifestly transverse currents — that 
have no bearing on gauge invariance — can always be 
added to any given prescription. From a phenomeno- 
logical point of view, therefore, it is unavoidable that 
one seeks a prescription that works best in reproducing 
the data (see, e.g., discussion in [23]). A number of the 
existing gauge-invariance preserving prescriptions have 
already been applied in this manner in pion photopro- 
duction [12-14] as well as in electroproduction [24, 25] 
reactions. 

However, most of the existing calculations based on 
phenomenological dynamical models are actually not 
gauge invariant. In fact, they revert to a variety of ad 
hoc recipes for the sole purpose of enforcing current con- 
servation, 

K Mii = o > (!) 

when the production current M M is on-shell, but not the 
gauge-invariance condition expressed by the generalized 
Ward-Takahashi (WT) identity [17, 18, 26] 

k^M* = -{F^Sp+kQ^ 1 + S^QfSj-klFvT) 

+ &;- P ' +k Q«\-P'\F t T) , (2) 

which is an off-shell condition. [This equation is repeated 
as Eq. (8) below, where its details are explained.] 

One of the few exceptions to this situation is the re- 
cent work by Pascalutsa and Tjon [14] where a fully gauge 
invariant pion photoproduction amplitude has been con- 
structed based on the Gross-Riska prescription [15]. This 
approach has been extended and applied as well to pion 
electroproduction [25]. The prescription of Ref. [15] also 
has been applied to the nucleon-nucleon bremsstrahlung 
reaction [27]. The Gross-Riska procedure relies on the 
observation that vertices and propagators always enter 
in the combination (vertex x propagator) in a given re- 
action amplitude. Therefore, a vertex form factor that 
depends only on the momentum of the propagating (off- 
shell) particle can be incorporated into the corresponding 
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propagator instead of being associated with the vertex. 
Of course, if more than one leg of the vertex belongs to an 
off-shell particle, this restricts the phenomenological form 
factors to being separable functions of the respective leg 
momenta. In addition, the form factors should be such 
that they do not lead to unphysical behavior of the re- 
sulting propagators. Gauge invariance is then fulfilled by 
constructing electromagnetic currents that obey the WT 
identities with the respective hadronic propagators mod- 
ified by the inclusion of the form factors as described. 1 
At the tree level, this prescription completely removes 
the form factors from the longitudinal part of the reac- 
tion amplitude; i.e., only manifestly transverse parts of 
the production current carry any form factor dependence. 

In the present work, we construct a photoproduction 
amplitude based on the field-theoretical approach given 
by Haberzettl [18]. The full formalism is gauge- invariant 
as a matter of course. However, in view of its complex- 
ity and high nonlinearity, its practical implementations 
require that some reaction mechanisms need to be trun- 
cated and/or replaced by phenomenological approxima- 
tions. Our objective here is to preserve full gauge invari- 
ance, in the sense of Eq. (2), for this approximate treat- 
ment, but allowing for the presence of explicit hadronic 
FSIs. This problem has been treated already in Ref. [28] 
as a two-step procedure where the gauge-invariant treat- 
ment of explicit FSIs was added on to an already gauge- 
invariant tree-level amplitude that had been constructed 
according to the prescriptions given in Ref. [19]. The 
present approach instead starts from the full amplitude 
and derives a single condition for the mechanisms to be 
approximated that follows directly from the generalized 
WT identity (2). It thus is more general and not tied to 
any particular tree-level treatment. Moreover, we present 
a general scheme that allows a systematic way of includ- 
ing more complex reaction mechanisms into the proce- 
dure. At the lowest order, it is found that the essential 
aspects of the results found in Ref. [28] remain true. In 
contrast to the prescription of Ref. [15], the present ap- 
proach does not impose any restriction on the type of the 
hadronic form factors that can be used. Furthermore, the 
longitudinal part of the resulting reaction amplitude re- 
tains these form factors even at the tree-level. 

The present paper is organized as follows. In Sec. II we 
present our approach to construct a fully gauge invariant 
photoproduction amplitude. In Sec. Ill we illustrate the 
approach developed in sect. II by applying it to the pion 
photoproduction. Section IV contains a summary with 
our conclusions. Some details of the present approach as 
well as of its model application are given in the appen- 
dices. 



Note that the electromagnetic vertices constructed in this way 
in Refs.[14, 27] differ from each other by a transverse piece. 



II. FORMALISM 

In the following, for definiteness, we will explicitly 
consider the production of pions off the nucleon, i.e., 
7 + N — > 7r + N, but the formalism will of course apply 
equally well to the photoproduction (or electroproduc- 
tion) of any pseudoscalar meson. Moreover, at interme- 
diate stages of the reaction, we will ignore other mesons 
or baryonic states since they are irrelevant for the prob- 
lem at hand, namely how to preserve gauge invariance in 
the presence of FSI. 

As mentioned, our approach is based on the field- 
theory formalism of Ref. [18]. For the present purpose, 
however, we do not need to recapitulate its full details. 
Instead, we employ the summarizing diagrams of Figs. 1, 
2, and 3. As is well-known [1], the photoproduction cur- 
rent M M can be broken down according to 

M " = + M£ + M t " + M£ t , (3) 

where the first three terms describe the coupling of the 
photon to external legs of the underlying 7riViV vertex 
(with subscripts s, u, and t referring to the appropri- 
ate Mandelstam variables of their respective intermedi- 
ate hadrons). These terms are relatively straightforward 
and easy to implement in a practical application. How- 
ever, the last term, the interaction current M^ t , where 
the photon couples inside the vertex, explicitly contains 
the hadronic FSI; its structure is, therefore, more com- 
plex than that of any of the first three terms. We read 
off the diagrams enclosed by the dashed box of Fig. 1 (a) 
that 

ML = < are + U^G \Fr) 

+ XGo (M£ + M? + m£ are + U"G Q \Ft)) , (4) 

where m£ are is the bare Kroll-Ruderman contact current, 
subsumes all possible exchange currents (see Fig. 2), 
Go describes the intermediate irN two-particle propaga- 
tion, and the FSI is mediated by the nonpolar part X 
of the ttN T matrix. Following Ref. [18], the notation 
\Ft) is used for the dressed N — > ttN vertex (including 
its full coupling-operator structure); for Ntt — ► N with 
the pion leg reversed, we use (Ft\. The isospin operator 
r (with its component index suppressed) is pulled out of 
the vertex explicitly for later convenience [see Eq. (12)]. 
The vertex obeys the equation 

|^r) = (l + XGo)|^barcr) , (5) 

where \F^ alc T) denotes the bare irNN vertex. 

Note that the ./V./V7 electromagnetic vertices appearing 
in Fig. 1(a), are also fully dressed vertices given by 

^ N = ^ Nhaic + m^ c G \FT) 

+ (Fr|G (m^ arc + M? + M£ + U»G \Fr)) , 

(6) 
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FIG. 1: (Color online) Diagrammatic summary of the field-theory formalism of Ref. [18]. Time proceeds from right to left, 
(a) Meson production current . The first line corresponds to Eq. (3) summing up, in that order, the s-, u-, and i-channel 
diagrams and the interaction current M^ t . (The different colors of the hadronic three-point vertices identify the wNN vertex 
in different kinematical situations.) The dynamical content of Mf^ t is explicitly shown by the diagrams enclosed in the dashed 
box of the last two lines. This also includes, in the bottom line, the final-state interaction mediated by the nonpolar 7riV 
amplitude X that satisfies the integral equation shown in (b). The diagram element labeled U subsumes all exchange currents 
contributing to the process (see Fig. 2). The diagram with open circle depicts the bare current fnJ^ (i.e., the Kroll- 
Ruderman term). (b) Pion-nucleon scattering with dressed hadrons. The full 7riV-amplitude is denoted by T, with X 
subsuming all of its nonpolar (i.e., non-s-channel) contributions. The latter satisfies the integral equation X — U + UGqX 
depicted in the third line here, where the driving term U sums up all nonpolar irreducible contributions to 7riV-scattering, i.e., 
all irreducible contributions which do not contain an s-channel pole (see Ref. [18] for full details). — In both parts (a) and (b), 
diagram elements with open, unlabeled circles describe bare quantities, and solid circles (or circles filled with colors) denote 
the corresponding dressed vertices and propagators. 
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FIG. 2: (Color online) Exchange-current contributions sub- 
sumed in ?7 M grouped by the topological properties of the 
underlying nN irreducible hadron contributions: The dia- 
grams subsumed in _B M are based on the photon attaching 
itself internally in all possible ways to a simple hadron ex- 
change graph, whereas D M subsumes the analogous contribu- 
tions arising from a triangular hadron graph; more complex 
structures are not shown explicitly. Note that the second and 
third diagrams of , and the first and third diagrams of 
explicitly contain the full off-shell interaction current M^ t . 
Implicitly it is contained in many more places. (The mid- 
dle diagram in the second row of _D M contains the diagram 
X^ where the photon is attached to the internal structure of 
the non-polar nN amplitude X. The details of this mecha- 
nism are irrelevant for the present considerations; they can 
be found in Ref. [181.) 




FIG. 3: Diagrammatic representation of the nucleon's elec- 
tromagnetic current obtained by attaching a photon to the 
second line of diagrams in Fig. 1(b). The terms correspond 
to those of Eq. (6). 



and illustrated diagrammatically in Fig. 3. 
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denote the bare vertices for 7iV — > N and ^ttN 



is the Kroll-Ruderman current m^ aie with 



the pion leg reversed). Also, the nucleon propagator, S, 
illustrated in the second row from the top in Fig. 1(b), is 
fully dressed according to 



(F haieT \G \FT) 



(7) 



such that the WT identity as expressed by Eq. (10a) be- 
low is satisfied. Share stands for the bare nucleon propa- 
gator. 

As alluded to in the Introduction, the production cur- 
rent M M is gauge invariant if its four-divergence satisfies 
the generalized WT identity [17, 18, 26] 

k^M» = -\F s T)S p+k QiS- 1 + S^QfSp^klFuT) 

+ &-- p > +k Q«\-?\FtT) , (8) 

where p and k are the four- moment a of the incoming 
nucleon and photon, respectively, and p' and q are the 
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four-momenta of the outgoing nucleon and pion, respec- 
tively, related by momentum conservation p' + q = p + k. 
S and A are the propagators of the nucleons and pions, 
respectively, with their subscripts denoting the available 
four-momentum for the corresponding hadron; Qi, Qf, 
and Qn are the initial and final nucleon and the pion 
charge operators, respectively. The index x at [F x t] la- 
bels the appropriate kinematic situation for ttNN vertex 
in the s-, u-, or i-channel diagrams of Fig. 1. This is an 
off-shell condition. In view of the inverse propagators 
appearing in each term here, if all external hadronic legs 
are on-shell, this reduces to 

fc M AF = (hadrons on-shell) , (9) 

which describes current conservation. 

Physically relevant, of course, is only current conserva- 
tion. However, the reason one must satisfy the off-shell 
condition (8) for gauge invariance to hold true is the re- 
quirement to have consistency across all elements of the 
underlying reaction dynamics. In view of the inherent 
nonlinearity of the process (due to the fact that the num- 
ber of pions is not conserved), the elements contributing 
to the full amplitude couple back into themselves 
nonlinearly [18]: For example, as can be seen from Fig. 2, 
the sum of exchange currents internally also contains 
the interaction current M(^ t in several places, with at 
least one hadron leg off-shell even if all external hadrons 
are taken on-shell. It is then found that it is not possible 
to achieve current conservation consistently unless the 
current satisfies the off-shell condition (8), which trans- 
lates into the condition (11) for the interaction current 
given below. 

The electromagnetic currents for the nucleons and the 
pions, T^r and T^, respectively, satisfy the WT identities 

fc/JW,p) = S^Qn - QnS; 1 , (10a) 
k^(q', q) = - Q.A- 1 , (10b) 

where the four-momentum relations p' = p + k and 
q' = q + k hold. It is therefore possible to replace the 
generalized WT identity (8) by the equivalent gauge- 
invariance condition 

K M L = -\F s T)Qi +Q f \F u T}e f + Q*\F t r) 

= -\F.)ei + \F u )e f + \F t )e^ , (11) 

where the operators 

&i = TQi, e/ = Q/r, and e n — Q v t (12) 

describe the respective hadronic charges in an appropri- 
ate isospin basis (component indices and summations are 
suppressed here), i.e, apart from some numerical factors, 
the e x are essentially given by the charges of the respec- 
tive particles. Charge conservation for the production 
process then simply reads 




FIG. 4: Diagrammatic representation of Eq. (14). 

In the following, it is more convenient to use the condition 
(11), instead of (8), together with (10). 

We emphasize here that if the single-particle electro- 
magnetic currents satisfy the WT identities (10) and if 
the four-divergence of the interaction current is given by 
(11), then the corresponding production amplitude M M 
will be gauge-invariant as a matter of course even if prop- 
agators and vertices have been subjected to approxima- 
tions. 



A. Preserving gauge invariance 

The preceding considerations are completely general. 
In practical applications, however, one will not be able to 
calculate all mechanisms that contribute to the full reac- 
tion dynamics and one must make some approximations. 
This is particularly true for the complex mechanisms that 
enter as depicted in Fig. 2. Approximations should 
preserve the gauge invariance of the amplitude. To see 
how this can be done, let us define 

M£ = m!t m + U»G \FT) , (14) 

which is shown in Fig. 4, and thus write 

= M£ + XG (M£ + M? + M») , (15) 

and recast the gauge-invariance condition (11) as a con- 
dition for M£. This immediately produces 

fc M M£ = (1 - UG )[- \F s ) ei + \F u )e f + \F t )e„] 

-UGofaiMS + M?)], (16) 

where 

(1 + AGo)- 1 = 1 -UGq (17) 

and 

U={l + XG y 1 X (18) 

were used, with U being the sum of all non-polar hadronic 
driving terms [cf. third line of diagrams in Fig. 1(b)]. In 
the last term of (16) obviously only the non-transverse 
parts of the u- and i-channel currents Mff and Mf will 
contribute. Denoting those respectively by m(J and raj 1 , 
i.e., 

fc„(M£-r<)=0 and fc M (M t " - m?) = , (19) 
we finally have 

fc M M£ = (1 - UG ) [ - \F s ) ei + \F u )e f + \F t )e n ] 



(13) 



-k p UG {m<t + m1) 



(20) 
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as the necessary condition that M£ must satisfy so that 
M? nt yields the gauge- invariance condition (11). Equa- 



tion (20) is exact 
to here. 



no approximation has been made up 



1. Approximating M£ 

The structure of the preceding condition suggests the 
following approximation strategy. The condition evi- 
dently is satisfied if we now approximate M£ by 

M£ = (1 - UG )M? - UG Q (m£ + m?) + , (21) 

where Mjf can be any contact current satisfying 

k^M* = -\F s ) ei + \F u )e f + \F t )e v (22) 

and T M is an undetermined transverse contact current 
that is unconstrained by the four-divergence (20). With 
the choice (21), the corresponding approximate M^ t is 
then easily found from (15) as 



+ 



XG Q [(M£ <) + (M? m?) + T"] . (23) 

In this scheme, therefore, the choice one makes for M[? 
(and T^) corresponds to an implicit approximation of the 
full dynamics contained in the right-hand side of Eq. (14). 
Moreover, beyond this actual choice, the only explicit ef- 
fect of the FSI X is from explicitly transverse loop con- 
tributions, which is precisely the same result that was 
found in Ref. [28]. Thus it follows that 



(24) 



and this approximate interaction current then obviously 
satisfies the gauge-invariance condition (11). 

Equations (21) and (23), together with the prescrip- 
tions for and T 11 as given in the following two sub- 
sections, are the main results of the present work. 

Note that the choice of T M , while it has no bearing on 
the gauge invariance itself, will have a direct effect on 
how, if at all, the FSI enters the approximate treatment. 
For example, putting for the moment 



T M = -UG [(MZ - m£) + (M t " - m: 



simply provides 



(25) 



(26) 



Therefore, this particular choice completely eliminates 
the explicit occurrence of the FSI and, for phenomcno- 
logical choices of M^, such as Eq. (35) below, this corre- 
sponds to the tree-level approximation where the full in- 
teraction current is replaced by a phenomenological con- 
tact current. We emphasize, however, that we do not ad- 
vocate actually using Eq. (25). This particular (extreme) 
choice merely illustrates that the undetermined may 



contain pieces that may be capable of partially undoing 
the explicit inclusion of the FSI in (23). In Sec. II A3 we 
introduce a phenomenological procedure for obtaining T M 
from the data. 

With the approximation (21) for M£, the dressed 
NNj vertex (6) becomes 



r^-r^ (baro) +mf: arc Go|^r} 

+ (F barc r|G (M^ + m£ + m?) 

+ (Ft | G [ (M£ - m{|) + (M t " - <) + T" 



(27) 



where Eq. (5) has been used. To the extend that the 
approximated M£ (21) fulfills the same condition (20) 
satisfied by the exact current M£, the above NNj vertex 
also satisfies the same WT identity (10a) obeyed by the 
exact dressed vertex given by Eq. (6). 



2. Choosing 

The phenomenological choice that we make here for 
M£ is a variant of the procedure proposed in Refs. [18, 19] 
that is more general than what was suggested in [28] . We 
parameterize the ttNN vertices by 



A 



A) 



m + m' 



fx i 



(28) 



where x — s, u, or t indicates the kinematic context, g n is 
the physical coupling constant, m and ml are the nucleon 
masses before and after the pion is emitted/absorbed and 
the parameter A allows for the mixing of pseudoscalar 
(PS: A = 1) and pseudovector (PV: A = 0) contributions. 
For simplicity, the functional dependence f x of the vertex 
(which depends on the squared four-momenta of all three 
legs) is chosen as common to both PS and PV, and it is 
normalized to unity if all vertex legs are on-shell. We 
define then an auxiliary current 



C = 



6 7i 



(2q - ky 



t-q 2 

(2 P ' - ky 



{ft ~ F) 



r, (U-F) 

u - p' z 

{2 P + ky 



s — p A 



-ifs-F), 



(29) 



where 



F = 1 - h (1 - S s f s ) (1 - 6 u f u ) (1 - S t ft) . (30) 

The factors S x are unity if the corresponding channel 
contributes to the reaction in question, and zero other- 
wise. In principle, the parameter h may be an arbitrary 
(complex) function, h — h(s,u,t), possibly subject to 
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crossing-symmetry constraints. 2 (However, in the appli- 
cation discussed in the next section, we simply take h 
as a fit constant. Note that h = corresponds to Ohta's 
choice [16].) With this choice for F, the auxiliary current 
C M is manifestly nonsingular, 



- e/ (2p'-fc)"Azi 
u — p' z 



s — p z 



l-h(l-S s f s )(l-S u f u ) 

l-h(l-S,f a )(l-S t ft) 
l-h(l-6 u f u )(l-S t f t ) 



(31) 



i.e., it is a contact current, and in view of charge con- 
servation, e,r + e/ — ej = 0, its four-divergence evaluates 
to 



kfj.C = e w f t + e f f u - afs 



(32) 



With the vertex parametrization (28), the gauge- 
invariance condition (22) may now be written explicitly 



as 



A + (l-A) — 



7JV + 771 



ml + m 



(33) 



or, equivalently, as 



A + (l-A)4-^ 

TO' + TO 



- {1 - X] ^T^( eJs - eff ^f • (34) 

where the respective terms in the braces differ by a man- 
ifestly transverse term. We can exploit this ambiguity 
and set 



A+(l-A)^ 



TO' + TO 



-(1-A) — 



ml + to 



e*ft - 0k p C 



(35) 



where all terms depending on the free parameter (3 sum 
up to a transverse piece Tg cx /3(7 M k p C p ~^C^). The pa- 
rameter (3 then allows us to mix between the pseudovec- 
tor content' found in the expressions within the braces 



on the right-hand sides of Eqs. (33) and (34), which cor- 
respond to (3 = and (3 = 1, respectively. Note that 
(3 = amounts to a 'more traditional' treatment of the 
Kroll-Ruderman term where the bare 757^^ coupling is 
simply dressed by the ^-channel form factor e n — > e v f t - 
For (3 = 1, by contrast, this dressing occurs via the linear 
combination — > ej/ s — e//„ of s- and u-channel form 
factors (which in general would be non-zero even for n° 
production). 

Obviously, the above choice of Mjl is not unique, for 
we can always add another transverse current to it. In 
this sense, the pieces proportional to the parameter (3 
in Eq. (35) is just a particular choice of the transverse 
current added to the contact current. We emphasize, 
however, that the transverse contact current in M(l must 
not be confused with the transverse contact current T p 
appearing in Eq. (23). Note, in particular, that Mjl (and 
any of its contributing pieces) does not appear inside the 
FSI loop integral, but does. 



3. The transverse contact current T M 

The most general structure of a transverse contact cur- 
rent in pion photoproduction 3 can be written as [4] 



where 



T\ = -a^K = - (7^ - fc") , 

TO TO 



(36) 



(37a) 



T2 = -3 [P» (2q ■ k - k 2 ) - (2q - kf P ■ k] , (37b) 

T 3 = ^2 (V? ■ k - M) > ( 37C ) 

= -L ( 7 ^p . k _ _ T f , (37d) 

TO 

with P^ 1 = (p + p')^/2. The operators constitute a 
complete set of manifestly transverse operators for real 
photons. The coefficients Aj should be free of any sin- 
gularities in order to ensure that T p is a genuine contact 
current. The simplest approximation one can make for 
these coefficients is to assume them to be of the form 



(38) 



with fco denoting here the photon energy and aj being 
dimensionless constants to be fixed by the data. Notice 
that the factor l/fc in this equation will be canceled by 



2 Regarding crossing symmetry, note that the form (30) for F ad- 3 For simplicity, we consider here only real photons, but these con- 
drcsscs the concerns raised in Ref. [21] regarding the original siderations can easily be extended to electroproduction processes 

choice for F made in [19]. as well. 
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the factor ko which appears once the matrix element of 
the transverse current T M is calculated. 

To explore the energy range of where the approxima- 
tion (38) may be expected to produce reasonable results, 
let us calculate the corresponding matrix element of T^, 



(39) 



where denotes the photon's polarization vector and up 
is the nucleon spinor, with three-momentum p, normal- 
ized as uplift = 1. (Note that up here does not contain the 
Pauli spinor, i.e., T is an operator in spin- 1/2 space.) We 
have, in the center-of-momentum frame of the system, 

f = F x a ■ e + iF 2 e • n 2 + F 3 a ■ k e ■ q + F 4 a ■ q e ■ q , (40) 

with n 2 = (k x q)/\k x q\, where the hats denote unit 
three-vectors (i.e., v = v/\v\ for an arbitrary vector v), 
and 



Fi=a.Q+[a 1 + a 2 — cos 6 ) — cosf 
m I m 



F 2 



i — cos 


\ \<t\ . 
— sin 


m 


/ m 






m j 


' m 



F 4 = n 2 ^ 



(41a) 
(41b) 
(41c) 
(41d) 



with cos# = k ■ q. The quantities c^, Si, and r/i are 
given explicitly in Appendix A in terms of the coeffi- 
cients a,j of Eq. (38). These results show explicitly that 
the constant approximation for the coefficients aj of (38) 
leads to a transverse contact current that accounts 
for parts of the partial wave contributions up to D waves 
in the final wN state 4 . Therefore, the constant approx- 
imation (38) should be suited for energies not too far 
from threshold. For higher energies, if higher partial- 
wave contributions should be needed, one might expand 
Aj in terms of Legendre polynomials of higher order and 
fit the corresponding coefficients. 

The isospin structure of the transverse contact current 
can be included explicitly by writing the coefficient Aj in 
Eq. (36) as 



3 

E 

i=l 



(42) 



in which case A -, A- 



approximated by Eq. (38). 



and At individually are to be 



See also Ref. [29], where the coefficients Fj in Eq. (40) are given 
in terms of the partial-wave matrix elements to any desired order 
of the expansion. 



B. Explicitly incorporating the dynamics of 
exchange currents 

The fitting procedure of T M discussed in the previous 
subsection provides an indirect phenomcnological means 
of accounting for the transverse parts of the exchange- 
current contributions of subsumed in Fig. 2 which 
are neglected when approximating M£ as shown in Fig. 4. 
Specifically, none of the transverse mechanisms subsumed 
under the currents E^ and , etc., as defined in Fig. 2 
explicitly enters the approximation procedure described 
so far. 

This can be done, however, in a systematic order-by- 
order manner. To this end, note that with 

yv. = _gM + D n + . (43) 
as defined in Fig. 2, we may write Eq. (14) as 

M£ = E»G \Ft) + M'f , (44) 

where 

M'» = m£ aie + D»G \FT) + ... . (45) 

We may now subject M'^ to the same approximation 
procedure employed previously for Mjf, now, however, 
explicitly taking into account the current E^ . 

One easily finds that both Eqs. (22) and (23) remain 
valid, but T M is now given by 



T n = ( E i* - e ^)Go\FT)+T'^ , 



(46) 



where e M is the non-transverse part of E^ (the same way 
is the non-transverse part of M£, for x = u,t), i.e., 



k^E» - e»)G \Fr) =0 . 
The remaining transverse current T ,Ai , 
k^ = , 



(47) 
(48) 



remains undetermined by this procedure. It plays the 
same role, obviously, at the present level that T p played 
at the previous level of approximation. We may, there- 
fore, either treat it completely phenomenologically in the 
manner of Eq. (36), or we may, in principle, treat its dy- 
namics explicitly by now incorporating the triangle-graph 
currents £> M of Fig. 2. This then leads to 



T / M = ^ M -d^GolFr) + T ,,fM 



(49) 



where d M is the non-transverse part of D 1 * and T"^ the 
remaining unspecified transverse contribution. 

In principle, one may in this manner include more and 
more complex dynamical mechanisms explicitly into the 
formalism in a step-by-step procedure. In practice, how- 
ever, even incorporating the first step, Eq. (46), explic- 
itly in the interaction current (23) is a highly non-trivial 
task since this involves a double-loop integral which is 
very costly to evaluate numerically. 
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Note, however, that as far as gauge invariance is con- 
cerned, one need not work with the full currents or 
D 1 *, etc. Providing everything else is done consistently in 
the manner outlined here, any approximation of, for ex- 
ample, that satisfies the gauge-invariance constraint 

k^(p\q\p,q)^Q' 7[ E(p',q' -k,p,q) 

+ Q' N E(p' -k,q',p,q) 

- E(p',q',p,q + k) Q n 

-E(p',q',p + k,q)Q N (50) 

will preserve gauge invariance as a matter of course [18]. 
Here, E describes the simple exchange graph obtained by 
stripping the photon off any of the three contributions to 
E^, and q,q' and p,p' are the incoming/outgoing pion 
and nucleon momenta, respectively, related by 

p + q + k=p' + q', (51) 

and Qtt, Q' v and Qn, Q'n the corresponding charge oper- 
ators. A similar equation can be written down for 
and, for that matter, for any topologically distinct con- 
tribution to U 11 . 



C. Unitarity 

The full field-theoretical photoproduction formalism as 
summarized in Fig. 1 is unitary as a matter of course 
(within the limits of the usual one-photon approxima- 
tion). It is a straightforward exercise to show that the 
discontinuity contributions for the corresponding unitar- 
ity relation arise from the intermediate nN propagator 
Go in the hadronic relation 

X = U + UG Q X (52) 

for the non-polar nN amplitude and from the dressed 
single- nucleon s-channel pole term \Ft)S(Ft\ that ap- 
pears in the full nN T-matrix, 

T= \Ft)S(Ft\ +X . (53) 

(The preceding two equations are depicted in the last 
two lines of diagrams in Fig. 1(b); see also Ref. [18].) In 
the one-photon limit, therefore, only the cut structure of 
the hadronic part of the photo-reaction is relevant and 
other cuts do not contribute to the unitarity relation. 
This applies, in particular, to the cuts that are contained 
in the contact current M£ shown in Fig. 4. Hence, any 
approximation of M£ will preserve the unitarity of the 
photoproduction amplitude. The approximation choice 
(21), therefore, does not violate unitarity. 

III. APPLICATION: PION 
PHOTOPRODUCTION 

In this section, we apply the approach developed in 
the preceding section to the photoproduction reaction 



7 + N — > 7T + N. At this stage, this is intended to be 
more of a feasibility study rather than a serious attempt 
at describing all features of the data. We will therefore 
make some simplifying assumptions along the way. 

In the present application, we restrict ourselves to pho- 
ton energies up to about 400 MeV. Therefore, in addition 
to the basic nucleons and pions discussed in the preced- 
ing section, our model also incorporates intermediate A's 
in the s- and w-channcls. Details of the dressing of the 
electromagnetic NA transition vertex in the present ap- 
proach is given in Appendix B. We also include the p, u, 
and oti meson exchanges in the t-channel. Note here that 
transition currents between different hadronic states are 
transverse individually and therefore play no role for the 
issue of gauge invariance. The details of the respective in- 
teractions are specified by the Lagrangian densities given 
in Appendix C. Form factors are attached to the hadronic 
vertices to account for the off-shellness of the respective 
intermediate hadrons. The details of these form factors 
are also found in Appendix C. 

For the ttN FSI, we employ the irN T-matrix devel- 
oped by the Julich group [30] which results from a dy- 
namical model based on a coupled-channels approach. 
Among other things, this interaction fits the irN phase- 
shifts and inelasticities below about 1.5 GcV. For larger 
energies it provides a background due to final state in- 
teractions which has to be supplemented by baryon res- 
onances. It should be noted that the Julich irN in- 
teraction is based on time-ordered perturbation theory 
(TOPT) [31]. Therefore, to be fully consistent with 
this interaction, one should also evaluate the amplitudes 
(x = s,u,t) within TOPT. In the present appli- 
cation, however, we have ignored this consistency re- 
quirement and evaluate these amplitudes following the 
Feynman prescription (which coincides with TOPT at 
the tree level). As a consequence, there is an ambiguity 
in defining the zeroth component of the four momen- 
tum of the intermediate state in the transverse ampli- 
tudes (Mg — m£) and (M t M — mf) appearing under the 
loop integral of the FSI contribution in Eq. (23). We 
follow the choice made (based on the gauge invariance 
consideration) in Ref. [11] for the zeroth component of 
the intermediate particle momentum in evaluating these 
amplitudes. 5 

In the present application, for simplicity, we ignore the 
dressing effects in the s-channel nucleon pole propagator. 
Moreover, we ignore the explicit dressing of the nucleon 
electromagnetic vertex as given by Eq. (27). Instead, we 
take the vertex given by Eq. (C4a) with physical coupling 
constants. The analogous approximation is also adopted 
for the NA<y vertex given by Eq. (C4g). Here, the cou- 
pling constants G\ and G2 are treated as free parameters 



5 Note that, in the present approach, the photoproduction ampli- 
tude is gauge invariant independent of this particular choice of 
defining the zeroth component of the intermediate particle mo- 
mentum. 
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TABLE I: Model parameters fitted to the reaction 7 A — + nN. 
The AA7 coupling constants Gi and G2 are constrained by 
the measured E2 to Ml ratio of Rem = -2.5% [32, 33]. The 
parameter /3 in the contact current is fixed to be zero from 
the outset. Moreover, we consider only pure pseudovector 
coupling and therefore always have A = 0. 



Gi 


G 2 


A (MeV) 


h 


3.84 


-1.94 


604 


0.01 



adjusted to reproduce the data. Such an approximation 
is not critical for the present purpose of illustrating how 
the method developed in the previous section works. Ob- 
viously, dressing of the electromagnetic vertices is more 
critical in electroproduction processes. 

In the following, we list the free parameters of the 
model in the present application. 

1) Gi and Gi: The two dressed electromagnetic cou- 
pling constants at the NAj vertex, as specified in 
Appendix C, are not independent of each other, for 
we impose the E2 to Ml ratio to be Rem — —2.5% 
as determined by the Mainz group [32, 33]. 

2) h, p, and A: The parameters h and f3 appear in 
the contact current in Eqs. (31) and (35), and the 
PS/PV mixing parameter A appears in Eq. (28). 
We take = and A = (pure PV) from the 
outset, so that h is the only parameter to be fixed 
in the contact current Mg. 

3) T^: The transverse contact current in Eq. (23) is 
found to be negligible in the present application, 
i.e., the parameters Oj in Eq. (38) are taken as Oj = 
0, (j — 1, 2, 3, 4). In other words, we found no need 
for such a current in reproducing the cross-section 
data. 

4) A: This regularization parameter is needed for the 
loop integral in the FSI contribution. We regularize 
this integral by introducing a cutoff function 



where q denotes the relative momentum of the in- 
termediate ttN state in the loop integral. Of course, 
this regulator may also be interpreted as the form 
factor which accounts for the off-shcllncss of the 
pion and nucleon in the loop integral. 

With the considerations mentioned above, we are 
left with only three independent free parameters in the 
present model. They are adjusted to reproduce the pion 
photoproduction cross-section data. The resulting pa- 
rameter values are given in Table I. Note that h is nearly 
zero, corresponding practically to the Ohta's choice [16]. 

Figure 5 shows the total cross-section result for the re- 
action 7 + p — » 7T° +p from the threshold up to T 7 w 400 




T Y (MeV) 

FIG. 5: (Color online) Results corresponding to the parameter 
set in Table I for the total cross section as a function of photon 
incident energy T 1 in the reaction 7 + p — > it + p. The 
dashed curve corresponds to the Born contribution and the 
dash-dotted curve to the FSI loop contribution. The solid 
curve is the total contribution. The data are from Ref. [34]. 



MeV. As we can see, the agreement with the data is very 
good except for energies above T 7 ~ 360 MeV, where the 
prediction tends to underestimate the data. In particu- 
lar, around T 7 w 390 MeV, the discrepancy is about 10%. 
We also see that the FSI loop contribution is relatively 
small compared to the Born contribution. However, it 



o + 

y+p— >ji +p y+p-->ju +n y+n— >it -+p 




n 1 1 1 - - 1 1 1 1 1 1 1 1 1 1 1 1 

60 120 60 120 60 120 180 

e (deg) 

FIG. 6: (Color online) Results corresponding to the parame- 
ter set in Table I for the differential cross sections in the cm. 
frame of the system in the reaction 7 + A^7r + Aat vari- 
ous photon incident energies T 7 . In the top row, the dashed 
curves correspond to the results represented by solid curves 
multiplied by an arbitrary factor of 1.1 The data are from 
Refs. [33, 35]. 
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plays a crucial role in reproducing the observed energy 
dependence through its interference with the dominant 
Born term. 

Figure 6 shows the results for differential cross sections 
for neutral and charged pion productions at various en- 
ergies together with the data. We see that, overall, the 
data are reproduced quite well. The dashed curves in 
the top row correspond to the results represented by the 
solid curves multiplied by an arbitray factor of 1.1. They 
are shown here to facilitate visualizing that the shape of 
the angular distribution is well reproduced, in spite of 
the absolute normalization being underestimated at this 
energy by ~ 10%, as can be seen better in Fig. 5. 



IV. SUMMARY 

By exploiting the generalized Ward-Takahashi identity 
for the production amplitude and total charge conserva- 
tion, we have constructed a fully gauge-invariant (pseu- 
doscalar) meson photoproduction amplitude which in- 
cludes the hadronic final-state interaction explicitly. The 
method is based on a field-theoretical approach devel- 
oped earlier by Haberzettl [18]. It is quite general and 
can be readily extended to any other meson photo- and 
electroproduction reactions. This method should be par- 
ticularly relevant for the latter reaction. 

As an example of application of the present approach, 
we have calculated both the neutral and charged pion 
photoproduction processes off nucleons up to about 400 
MeV photon incident energy which illustrates the feasi- 
bility of the present method. 

Obviously, for a more quantitative calculation, includ- 
ing not only cross sections but also other observables, 
some of the approximations made in the present fea- 
sibility study should cither be improved or altogether 
avoided. In particular, the dressing of the electromag- 
netic vertices as given by Eq. (6) needs to be carried out. 
Also, for pion photoproduction, one should constrain the 
parameters of the present model, if possible at all, by 
comparing the amplitude of the present approach in the 
chiral limit with that of the Chiral Perturbation Theory 
Work in this direction is in progress. 



are 



a = N 



aiG + a 3 



m \ m 



(Ala) 



ai = N [- ai JG - a 3 (1 + uj q H) + a 4 (1 - IH - 2JG)} , 

(Alb) 

a 2 = N [a 3 — a 4 ] mH , (Ale) 



where u> q = \J q 2 + m 2 and 
N E 
Gee 1 



e q + m jek+m 1 
2m 

1*1 



2m m 



H 



Sk+m 
\k\ 



(e q + m) (ek + m) 
/ eee e k + \k\ +e q , 



J = 



m 



(A2a) 
(A2b) 

(A2c) 

(A2d) 
(A2e) 



e q + m 

with e p ee yjp 2 + m 2 . Furthermore, 

0i=N [-ai JG + a 3 uj q H + a 4 {IH + 2JG)] , (A3a) 
02 = -a 2 , (A3b) 

and 



5i = N 



ai JG-a 2 ^K + a 3 (G + Lu q H) (A4a) 



+ a A {G-IH- 2JG) 
t?2 = — N [a 3 + 04] mH , 
with K = 2{sk + m)/{s q + m), and finally 



m = N 



C12K — \ a 3 — 04 — (J — mH) 
\ mi 



(A4b) 



(A5) 



APPENDIX B: TRANSVERSALITY OF THE 
DRESSED iVA 7 VERTEX 
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APPENDIX A: TRANSVERSE CONTACT 
CURRENT 

In this appendix, we give the explicit formulas for the 
coefficients ai, 0i, Si, and m appearing in Eq. (41). They 



For completeness, in this appendix we show how to 
dress the N A7 transition vertex in the present approach 
and demonstrate that the resulting vertex is purely trans- 
verse as it should be. 

Bare transition currents can easily be made transverse 
by expanding the current in an appropriate transverse op- 
erator base. It is not obvious, however, that the transver- 
sality will remain true after one dresses the current. It 
will be shown here that this is indeed the case. As a spe- 
cific example, we will treat the electromagnetic current 
for the ./V — > A transition. 

The dressing mechanism for this current is depicted 
in Fig. 7 which is constructed in analogy to the dressed 
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nucleon current in (6) (see also Ref. [18] for full detail). 
We will consider here only dressing mechanisms in terms 
of nucleons and pions. Other particles are independent 
from the ones considered here and formally would not 
add anything new except complicating the presentation. 
The two equivalent forms arise from attaching the photon 
in all possible ways to two equivalent bubbles shown in 
Fig. 8. It should be obvious, however, that as a physical 
process the transition N — > A as shown in this figure is 
not possible because of isospin conservation. Within the 
present context, therefore, the bubbles in Fig. 8 form but 
the topological backdrop against which the current shown 
in Fig. 7 is constructed but are not considered as having 
a physical meaning of their own. 

Following Ref. [36] , the most general transverse current 
for 



equation 



7 (fc) + N(p) -> A(p') 



(Bl) 



may be written as 



= G x 75 (k p Y - gP" If) 

+ G 2 y 5 (kf } P»-g^k-P) 
+ G 3l5 (k> 3 k»-g^k 2 ) 



(B2) 



where P = (p + p')/2 = (2p + k)/2; the Lorentz indices /i 
and pertain to the photon and the A, respectively. The 
Gi are the corresponding form factors; for real photons, 
only G\ and G*2 contribute. This is the ansatz that one 
chooses for the bare current. 

To show that the transversality of the bare current is 
preserved if one now dresses the current, we use the first 
form of the current shown in Fig. 7. Using the nota- 
tion of Ref. [18], we can translate this into the schematic 




FIG. 7: (Color online) Dressed transition current for *yN — > 
A. Solid circles depict dressed vertices and currents, whereas 
open circles show the corresponding bare quantities. The two 
equivalent forms differ by whether the initial irNN vertex or 
the final tvNA vertex is fully dressed. Contributions with 
intermediate particles other than pions and nucleons are not 
shown. 



(f»\S N oA % \F) 
+ (f\S N o A^M^) 
+ {f\S N o[A 7r T^]\F) 
+ (f\[S N T%S N ]oA v \F) 



(B3) 



where the order of terms is exactly as in Fig. 7. The 
(transverse) bare current is denoted by T^ 1 ; f 13 is the 
bare irAN vertex and the corresponding bare contact 
current. The latter are given by 



f = fA 0^ qv , 



(B4) 



where q v is the incoming pion momentum, B^ v the cou- 
pling tensor, and /a the bare coupling constant, and by 



(B5) 



where e T is the charge of the intermediate pion. F is the 
dressed irNN vertex and M(^ t the corresponding dressed 
interaction current. The convolution of the intermedi- 
ate pion and nucleon propagators (with and without at- 
tached photons) is denoted by A o B. The momentum 
dependence is suppressed here, but can easily be found 
by noting that the initial nucleon momentum is p and 
that the photon feeds a momentum k into the line (or 
vertex) to which it is attached. 

Using the Ward-Takahashi identities for the pion, the 
nucleon, and the interaction current, 

k^SNT^SN = e' N (•S N - S N ») , (B6a) 

fc M A T r£A w = (.A, - A ff «) , (B6b) 

K M L = -F[»e N ] + [e' N »]F + [e ff .]F , (B6c) 

respectively, where the solid dot (•) indicates at which 
point the photon momentum is injected into the equa- 
tions. For example, 

e' N (.S N - S N *) = e' N \s N (p") - S N (p" + k)] , (B7) 



FIG. 8: (Color online) Topological structure of the bubbles 
underlying the construction of the transition current in Fig. 7. 
The first equality follows from the fact that the dressing for 
both bubbles is done in terms of the non-polar nN T-matrix 
X and both bubbles with dressed vertices, therefore, can be 
written in terms of the graphs of the second line. Taken as a 
physical process, the transition N — » A is not possible either 
on or off-shell because of isospin conservation. 
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where p" is the initial momentum of the nucleon with 
charge e' N within the loop. The notation [e^*] specifies 
that the photon momentum is injected into the (incoming 
or outgoing) particle line with charge e x . We now find 

k^* = k^ -fAeA0^^\S N o A„\F) 



-e N (f\S N oA n \F). 

+ e w (f \S N °[A n .]\F) 

+ e' N (f\[S N .}oA 7r \F) 

+ e jr (f\S N o[.A v ]\F) 

-e^(f\S N o[A^]\F) 

+ e' N (f\[,S N ]oA 7I \F) 

-e' N {fe\[S N *] A*\F) 

= -f^{e^k„\S N oA n \F) 

-e N (f\S N oA„\F)* 

+ e„(f\S N o[*A v ]\F) 

+ e' N (f \[*S N ]oA 7V \F) , (B8) 

In the second term, the dot • simply indicates that the 
overall four-momentum of (/^I-Sat A^F) is p+k. In the 
first, third, and fourth terms, the intermediate pion prop- 
agator depends on the same loop variable q in all terms, 
but the left-most vertices have momentum dependencies 
that can be combined according to 

-fAe v (O ^\S N o An\F) + e v (f\S N o [.A W ]|F> 

+ e' N (f\[.S N ]oA^\F) 

= fA[-e v (e^\+e v {e^(q + kU 

+ e' N {e^q v \]s N oA„\F) 

= fA[(e' N +e„){9^q v \\S N oA v \F) 
= (e' N + e n )(f\S N oA n \F) 

= e N (f\S N o A n \F) , (B9) 

where charge conservation, 

ejv = e T + e' N , (BIO) 

was used. Comparison with the first bubble of Fig. 8 
shows that 



Zan = (f\S N o A V \F) 



(Bll) 



is just equal to this topological bubble. Hence we have 



k^ 



e N 



j an 



(P) 



J AJV 



(p + k) 



As explained above, the S^jv do n °t describe a physical 
process and vanish individually. Hence, we have 



k^T^ = 



(B13) 



and the dressed current thus is also transverse. 

We observe that the present approach for dressing the 
NAj vertex differs from the dressing mechanism em- 
ployed in other approaches (see, e.g., [13, 14]) by the 
presence of the second diagram on the right-hand-side of 
the equality in Fig. 7. Note that this term involves the 
three-particle to one-particle transition, jttN — > A, and 
therefore is outside the model space considered in those 
approaches. We emphasize, however, that the presence 
of the AttNj contact vertex is absolutely necessary for 
preserving the gauge invariance of the dressed vertex in 
view of the momentum dependence of the NAir vertex as 
given by Eq. (C2e) below. This suggests, of course, that 
indiscriminate truncation of the model space along par- 
ticle numbers is not a good dynamical ordering scheme 
as far as the gauge-invariance condition is concerned. 



APPENDIX C: INTERACTIONS 

Our model for the s-, u-, and t-channel amplitudes 
Mf, M£ and Mf, respectively, in Eq. (3) is constructed 
from the interaction Lagrangian density written as a sum 
of two terms, £; nt = £hadr + £eicc, where £hadr denotes 
the part of the interaction Lagrangian containing only 
the hadron fields, and £ icc contains the electromagnetic 
interaction with hadrons. For £hadr, we have 

Aiadr = C-N Ntt + C*NN P + ^NN ai + £-NNu + £>N Att , (CI) 

with 

1 - A 



£-NNtv = -g-n^ 75 



Cnn p = * ( 7m - -P-°i»>V ) r- r * , (C2b) 



m + ml 



t? • f * , (C2a) 



2m n 



Cnnu = — * ( 7m 



C-NNcn 
£-NAtt 



2lTlN 



g ai *7a*75T- * , 

ft • (0„tF) * + H. c. 



(C2c) 
(C2d) 
(C2e) 



(B12) 



Here, \& and 'J'^ denote the nucleon and A fields, respec- 
tively, 7r the pion field, the p-meson, w M the w-meson, 
and Si the ai-meson fields. The latter is included as the 
chiral partner of the p-meson. The vector notation refers 
to the isospin space. stands for the isospin 1/2 to 3/2 
transition operator. The nucleon and pion masses are de- 
noted by rriN and m v , respectively. The g„, g p (K p ), g ai , 
9u>(kuj) and /nAtt are the corresponding coupling con- 
stants. We use g-x = 14.4 and JnAtt — 0.36 [30]. For 
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the NNp coupling constants, we use (g p /2) 2 /4w = 0.91 
and k p = 6.1 [37], while {g^/2) 2 /4tt = 11 and n u = 
[38]. The coupling constant g ai — m^f^/m^, with 
m ai w 1260 MeV denoting the mass of the ai-meson, 
has been fixed from the chiral-symmetry considerations 
following the work of Wess and Zumino [39] . 

The electromagnetic interaction Lagrangian density is 
given by 



£clcc 



(C3) 



with 



C NNl = -e * ( e 7p - -^o»vd v ) A" * (C4a) 



f 

Cnnkj = e — *7 57AI [f x tt] 3 VA^ , 
£™ 7 = e P^tt) x 7t] 3 ^ , 

777.77- 



(C4b) 

(C4c) 
(C4d) 

(C4e) 



x ^2[(0"tt) x a\ - {d"w) x 3$*] + tt x af^ , 

(C4f) 



+ e^-(^)Tt 75 ^+ H.c. , 



(C4g) 



where F pv = d^Ay — d u A p with A^ denoting the elec- 
tromagnetic field and atf" = d^a\ — d^a^. e is the pro- 
ton charge; e = (1 + t z )/2 and k = [1.79(1 + t z )/2 - 
1.93(1 — t z )/2] are the charge and magnetic moment op- 
erators of the nucleon, respectively. The coupling con- 
stants gw 7 = 0.374 and g plxl — 0.125 are fixed from the 
decay of the w- and p-meson into ir°+j, respectively. The 
signs of these coupling constants are consistent with those 
determined from the study of pion photo-production in 
the 1 GeV region [40] . e pa \ u is the totally antisymmet- 
ric Levi-Civita tensor with g 0123 = +1. The Lagrangian 
£ ai Tvy is obtained from £ ai7rp in Ref. [39] by combining 
it with the vector-dominance model. 



and are 



The propagators required for constructing Mjf, M£, 

(C5a) 



S n (p) 



1 



q 2 - m 2 . 



g» v - q^q v /m 2 v 



q 2 - ml 



for v = p, u), ai , 



1 



p 1 - m N 

i> + m A 
p 2 — m\ 



(C5b) 
(C5c) 



1 2p^p u p^Y - p v Y~ 

y 3 ' ' 3 m\ 3m A 

(C5d) 



where A w (q) denotes the pion propagator and D% v (q) the 
vector {p, u>) and axial-vector (ai) meson propagators. 
Sn(j>) and S^{p) are the nucleon and Rarita-Schwinger 
A propagators, respectively; m& — 1232 MeV denotes 
the mass of the A. Note that for the s-channel A reso- 
nance contribution, we have used the dressed VA7T ver- 
tex and the dressed A propagator according to Ref. [30] 
and consistent with the ttN FSI used. 

The amplitudes MJ*, M£, and constructed from 
the preceding Lagrangians are diagrammatically repre- 
sented in Fig. 1(a). 

Our model for M[f, Mj^ and M t M is supplemented with 
hadronic form factors, except for the s-channel A contri- 
bution where the dressed vertex is used. So, the NNir 
vertex in the s- and u-channels and the NAn vertex in 
the u-channel are multiplied by a form factor 



F B {p 



A% + (p 2 + m%f 



(C6) 



where p denotes the three-momentum of the off-shell 
baryon. In the above equation B stands for either the nu- 
cleon or A in the intermediate state. We take Kb = 1.2 
GeV for both baryons. 

The hadronic vertices in the t-channel Mj* amplitude 
are also supplemented by form factors of the form 



F a (q 2 ) = 



A 2 



A 2 



(C7) 



where a — tt, p, ui, a\. We take A T = 900 MeV and n T 
I [38] and A Q = 1850 MeV and n a = 2 as a — p,u),a\ 
[37]. 
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